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Abstract 

It has previously been shown theoreticahy that the exchange of hnear momentum between the 
hght field in an optical cavity and a vibrating end mirror can entangle the electromagnetic field 
with the vibrational motion of that mirror. In this paper we consider the rotational analog of this 
situation and show that radiation torque can similarly entangle a Laguerre-Gaussian cavity mode 
with a rotating end mirror. We examine the mirror-field entanglement as a function of ambient 
temperature, radiation detuning and orbital angular momentum carried by the cavity mode. 

PACS numbers: 03.65.Ud, 42.50.Pq, 42.50.Lc, 45.20.dc 



1 



I. INTRODUCTION 



The optomechanical driving of optical cavities with laser radiation is an experimentally 
promising method for eliciting quantum mechanical behavior from classical objects 
y, 0] . In the case of a two-mirror cavity with a perfectly reflecting end-mirror mounted on a 
cantilever and allowed to vibrate along the cavity axis, laser radiation entering through the 
slightly transmissive fixed input mirror can both trap and cool the mirror motion, depending 
on whether it is detuned slightly above or below a cavity resonance j^. Radiation pressure 
can thus decrease the number of phonons in the vibrating mirror. 



fijjj. 



eff 
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since cooling lowers the effective temperature of the mirror Tpg while trapping increases its 
effective vibration frequency cjcfr- A key goal is to achieve rim < 1, i.e. to place the mirror 
near its quantum mechanical ground state. Evidence of reaching that regime could come 
from measuring ra^, and a strategy for doing so has recently been proposed 6|. 

Radiation pressure can play an additional important role as a mechanism for generating 
entanglement. In addition to its basic interest, quantum entanglement - the existence of 
correlations disallowed by classical physics between two objects - has been demonstrated to 
be a valuable resource for quantum information processing tasks such as teleportation and 
superdense coding (see Ref. j3] and references therein). 

It has been shown theoretically that radiation pressure can entangle a light field to a 



vibrating mirror both in the absence [8| and in the presence 



entangle two vibrating mirrors, both when they form a cavity; 
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131 ]. A single vibrating mirror can also entangle two 



911 of a cavity. It can also 
~l 

111 ] and when they do not 



151 ] or more 16| light fields. The 



teleportation of quantum states from a light field to a vibrating mirror has been proposed 
in Ref. [17]. 

In these proposals the physical mechanism behind the generation of entanglement is 
the exchange of linear momentum between the photons in the cavity and the mirror. In 
recent work, we have proposed a rotational analog of mirror cooling and trapping where the 
light-mirror interaction relies on the exchange of angular momentum between a Laguerre- 
Gaussian cavity mode and a spiral phase element used as a cavity end-mirror 18|]. This 
article further investigates that system and shows how the radiation torque can generate 
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FIG. 1: (Color online). Arrangement for entangling the rotation of a mirror with a Laguerre- 
Gaussian (LG) mode using angular momentum exchange. A Gaussian (G) field of optical charge 
enters the cavity through a fixed partially transparent input coupler (IC) that does not affect its 
charge. A perfectly reflecting rear mirror (RM) rotating about the cavity axis on a support S adds 
a charge 21 to the beam on reflection. Both the IC and the RM are spiral phase elements. The 
angular deflection of the RM from equilibrium (<^o = 0) is indicated by the angle <j). The charge 
on the beams at various points is also indicated. 

entanglement between the optically charged light field and a rotating mirror. We investigate 
the entanglement as a function of ambient temperature, detuning of the laser radiation, and 
orbital angular momentum carried by the cavity mode. 

The rest of the article is organized as follows. Section [III briefly reviews the cavity config- 



uration proposed in Ref. ISj. Section [XTTl discusses the steady state of the system, discusses 
the linear fluctuations of the field and mirror motion about that steady state, and examines 
the dynamical stability of the rotating mirror. Section [IV] demonstrates the onset of entan- 
glement between the intra-cavity Laguerre- Gaussian mode and the rotating end-mirror and 
discusses its dependence on ambient temperature, detuning, and angular momentum of the 
light field. We also briefly discuss the experimental measurement of the entanglement. Sec- 
tion |V] provides a summary and a conclusion. Appendix Rl derives and solves the Lyapunov 
equation for the correlation matrix determining the entanglement between the mirror and 
the cavity mode. 

II. CAVITY MODEL 

The optomechanical system under consideration is described in detail in Ref. [l^ and 
illustrated in Fig. [H Both cavity mirrors are spiral phase elements [l^ . The input mirror 
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is slightly transmissive but does not change the orbital angular momentum of any beam 
passing through it. In reflection, however, it removes a charge 21 from the beam. Likewise 
the end mirror, which is assumed to be perfectly reflecting, is designed to add a charge 21 to 
a beam on reflection. With these specifications it can be seen that a Gaussian beam (with 
no charge) incident on the cavity can transfer a torque = clh/L per photon to the end 
mirror, where c is the velocity of light and L is the length of the cavity 18| . 
In dimensionless units, the Hamiltonian that describes this system is 18 



H = hiUca^a + ^{Ll + 0^) - hga^a<p. (2) 

where a and are the bosonic annihilation and creation operators for the cavity mode of 
frequency Uc, Lz and are the angular momentum and angular displacement, respectively, 
of the rear mirror about the cavity axis, with [Lz, (p] = —i. Finally a;^ is the angular rotation 
frequency of the rear mirror and 

is the optorotational coupling parameter, with / = MR'^/2 the moment of inertia of the 
mirror of mass M and radius R about the cavity axis. 



III. QUANTUM DYNAMICS 

The evolution of the operators describing the dynamics of the intracavity field and the 
mirror motion is given by quantum Langevin equations accounting for noise and damping 
due to the vacuum fluctuations entering the cavity field as well as to the Brownian noise 
coupling to the moving mirror [l^, 

a = — g(j))a — — a + 1/70™, 

= UJ^Lz, 

Lz = -u^<P + ga^a - ^Lz + (4) 

where 5 = — is the detuning of the laser frequency ul from the cavity resonance, 7 
is the cavity damping rate and is the intrinsic damping constant of the rotating mirror. 
The noise operator a" describes the laser field incident on the cavity, of mean amplitude 
(a™(t)) = a™. Its fluctuations are taken to be delta-correlated, 

{6a'''{t)6a''''\t')) = 6{t - t') 
4 



and add vacuum noise to the cavity modes. The Brownian noise operator e™ describes the 
mechanical noise that couples to the mirror from the environment. Its mean value is zero 



and its fluctuations are correlated at temperature T as [20| 



(5) 



£±- r ^e-(*-')^ [l + coth ( 

where k-Q is Boltzmann's constant. For high mechanical quality, /D^ ^1), the fluctu- 
ations become delta-correlated and Eq. ([5]) simplifles to 

{5i\t)5i-{t')) = ^{2n + mt-t'), (6) 

where 

is the mean number of thermal phonons available at the mirror of effective rotation frequency 
Ci;efT and temperature T [see Eq. ( I29l) for an expresion for cUefr]. The thermal phonon spectrum 
peaks approximately at the temperature 

We show later on that it is possible to maintain a mirror-fleld entanglement close to its 
maximum value for T < Tc, but that that entanglement decreases sharply beyond T T^. 
(Fig. [3]). This can be understood intuitively by noting that for T > we obtain uj > ujcs 
and the mirror can be rotationally excited by thermal phonons, but this ceases to be the 
case for T < Tc , in which case u! < Ues- 

A. Steady state 

The semiclassical, steady-state solution of Eq.(jl]) is 

_ ^a^^[^/2-z{6-g<Ps)] 



q\a 

- — 



(7/2)2 + _ g^^Y 

2 



L,,s = 0. (9) 
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The steady state field amplitude is complex in general, but in simple situations where 
the incident field is kept constant throughout, it is possible to chose its phase such that ag 
is real, 



(10) 



W + - 9<PsY]'"' 

We restrict our discussion to that situation here. The steady-state is bistable for high 
enough incident field 2l|, |22|], but we assume in the following that an electronic feedback 



loop maintains the length of the cavity and hence removes bistability. This is a standard 



practice in current experiments on optomechanical cooling 
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B. Fluctuations 



To account for the quantum fiuctuations about the semiclassical steady-state we expand 
the Heisenberg operator a about the state in the usual way as a = + 5a, and similarly 
for the other operators. Linearizing the quantum Langevin equations of motion in the 
fiuctuations gives then 

5'X = -^6X + A6Y + v^<5X;^ 

6Y = -^6Y - AdX + Gd(j) + y^dY^"", 

6(f) = uj^6L^, 

SL, = + GSX - ^SL, + 5e"\ (11) 

where A = 6 — g(j)s is the effective cavity detuning, G = gas\/2 is the effective optorotational 
parameter, and we have redefined the field fiuctuations in terms of their quadratures as 

6Xa = {6a + 5a^)/V2, 

6Ya = {6a - 6a^)/iV2. 
Eq. (ITT]) can be recast compactly as 

u{t) = Bu{t) + n{t), (12) 
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where u(t) = {S(f>,SLz,SXa,SYa) is the vector of fluctuations, n(t) 
(0, 5e™, y/^6X^J^, y/jSY^'^) is the input noise vector, and 

/ \ 

-uj^ -D^/I G 

-7/2 A 

\G -A -7/2/ 

This matrix determines the dynamic stabihty of the physical system, and also a measure 
of the entanglement between its two subsystems, the intracavity field and the rotating mirror. 



B 



(13) 



C. Dynamic stability 

According to the Routh-Hurwitz criterion, the stability of the steady-state solution is 
assured if none of the eigenvalues of the matrix B has a positive real part. This requirement 
can be reframed as a series of inequalities that have to be obeyed by the matrix elements of 
B 23|. These inequalities, given by 



+ (7^ + 4A2)2)] + IQG^Pcu^^^A + WK^' + 47A2) 
+ 4Djj(4cjj72 + 7^ + 4G^u^A + 4j^A^) > 0, 
uj^{-f^ + 4A2) - AG^A > 0. 

have been numerically verified to hold for the parameters of this article. 

IV. ENTANGLEMENT 

We quantify the entanglement between the cavity mode and the rotating end mirror in 



(14) 



terms o 
Ref. 



the logarithmic negativity En, a measure of entanglement that was proposed in 



241 ] to descri 



of examples 



3e continuous Gaussian variables and has been evaluated in a number 
261 . I27I I. The logarithmic negativity can be expressed in terms of the 



B y, y, H- 



elements of the correlation matrix 



Cij = [{ui{oo)uj{oo) + Uj{oo)ui{oo))]/2. 



(15) 
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For the problem at hand, C is a 4 x 4 matrix that can be cast in the form 

R F 

where R, F and 5" are 2x2 matrices. In terms of these matrices the logarithmic negativity 



C 



(16) 



IS 



24| 



En = max[0, -ln(2r/-)], 



(17) 



where 



= 2~'/^a - [a' - 4\C\]'/'f\ (18) 

and a = li^l + IS*!— 2|F|, \R\ being the determinant of the matrix R. Quantum entanglement 
occurs for E^ > 0, i.e. for r]^ < 1/2. 

The bipartite entanglement between the cavity mode and the rotating mirror can be 
evaluated directly from the knowledge of the matrix B 9|. Integrating Eq. (fT2|) formally we 
have 

u{t) = M{t)u{0) + I dsM{s)n{t-s), (19) 
Jo 

where M(s) = exp(i?s). When the inequalities in Eq. (fT4|) are obeyed the system is stable, 
as we have seen, hence M(oo) = so that 



t^oo 



u{oo) = lim / ds M{s)n{t — s), 



(20) 



and the steady state is independent of the initial conditions m(0) as it should be. 

The stationary quantum fluctuations are zero-mean Gaussian processes fully character- 
ized by the 4x4 correlation matrix Cij which reads, with Eq. (1201) . 



/*oo /*oo 

a, = Y, ds ds'M,k{s)M,i{s')<^u{s - s') 
kJ Jo Jo 



(21) 



where 



^ki{s - s') = {{nk{s)ni{s') + ni{s)nk{s')))/2 

characterizes the stationary noise correlations. In the limit of a large mechanical quality 
factor, Eq. (jSD yields readily ^ki{s — s') = DkiS{s — s'), where 

/O \ 

D^{2n + 1)/I 

7/2 

\0 7/2/ 



D 



(22) 
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FIG. 2: Logarithmic negativity -Eat as a function of ambient temperature, for the parameters of 
Table H 



and Eq. ( !2Tll simplifies to 



C= I dsM{s)DM{sf, 
'o 



which can be integrated by parts to give 

BC + CB^ 



-D, 



(23) 



(24) 



see Appendix lAl Using Eqs. ( fT3l) and (l22l) . Eq. (IMl) can be solved analytically to yield the 
correlation matrix C, see Appendix [Aj The full solution is quite complicated and its form 
is not particularly instructive, hence it will not be reproduced here. 

We have evaluated numerically the logarithmic negativity Ej^j for various parameters. 
Our main results are summarized in Figs. [2}|4] for the parameter values of Table. [U where 
the mechanical quality is 

Qt> = (25) 



and F is the optical finesse. Our parameters are somewhat similar to the case of the linearly 
vibrating mirror of Ref. 9|, hence we make comparisons between specific figures in that 
article and the present one when appropriate. 



A. High temperature regime 

Fig. [2] shows En as a function of the temperature T of the mirror environment. Our 
results, which resemble those reported in Fig. [2] of Ref . [9|, indicate that quantum entangle- 
ment persists until roughly lOOK, that is, it should therefore be measurable in experiments 
in cryogenic environments. 

9 




FIG. 3: (Color online) Solid line: Deviation /S.E^ of from its maximum value Eq ~ 0.5 as 
a function of T. The graph shows the low temperature end of the curve. Dashed line: mean 
number of thermal phonons n [Eq.([7|)] at frequency oocs as a function of T. The vertical dotted 
line corresponds to the temperature Tc at which thermal phonons of energy larger than hiOcS start 
exciting the rotor, degrading the entanglement. Parameters of Table HI 

We remark that the temperature in Fig. [2]T is not the temperature of the mirror; as is well 
known the presence of red-detuned laser light in the cavity can increase the mirror damping 
from D(f) to D^s and thus lower the effective temperature of the mirror T^s significantly 
below T. More precisely, 



-efr 



T. (26) 



For the case of T = 50K, for example, we find Teg = lOyuK for the parameters of Table | 



where is (see Eq. (7) of Ref. [18|) 



LUc + (7/2)2 

[(7/2)2 + _ A)2] [(^/2)2 + (^ + A)2] • (27) 
Here u is the response frequency of the system, taken in our case to be Ues [see Eq.( l29l) ]. 



B. Low temperature regime 

Figure [3] is a blown-up view of the portion of Fig. [2] for low ambient temperatures. The 
entanglement is maximum at T ~ 0, as would be intuitively expected. It remains at that 
value as the temperature is increased until a temperature Tc where it undergoes a sharp 
decrease. At that temperature the peak frequency of the thermal phonons becomes larger 



10 



0.25 
0.2 
0.15 

En 

0.1 
0.05 

1 2 3 4 5 

A 

FIG. 4: The logarithmic negativity plotted as a function of the detuning A (normalized by the 
angular frequency uj^ of the rotating mirror) of the laser radiation from the cavity resonance. The 
three different curves correspond to different masses of the rotating mirror. The maximum of the 
entanglement occurs practically at A = i^cff ^ O.Sw^ (see the discussion in Section IIV Cp for all 
three curves. For this plot we used the parameters from Table |T] except for the following changes : 
L = 0.1mm, R = 250/im, = W, F = 10^, and T = 400mK. 

than tUcfT and they are capable of randomly exciting the rotating mirror, thus degrading its 
entanglement with the cavity mode. 

This point is further illustrated in Fig. [3l which shows the mean number of thermal 
phonons n at frequency Ucs, [Eq.([7])] as a function of T. It increases noticeably for temper- 
atures above and is clearly anti-correlated to E^. This is consistent with an approximate 
calculation that shows that to lowest order in n we have 

En ^Eo- nn, (28) 

where Eq is the maximum value of the entanglement and k is a constant whose value is a 
function of the parameters of the system. The analytical forms of Eq and k are involved. 
For the parameters of Table [T] their values are Eq ~ 0.5 and k ^ 3.1 x 10~^. 




C. Detuning 

Fig. m shows the logarithmic negativity as a function of the detuning A = ujc — ujl — 4^3 = 
for three masses of the rotafng mhxor. Th.s result ,„ay be eompared 
to Fig. 1 of Ref . [9( . A careful examination of the detuning results reveals that for each mass 
the entanglement is maximized when the detuning is equal to the mirror's effective angular 
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frequency of rotation A = Ucs, As is known, the presence of red-detuned radiation (A > 0) 
causes anti-trapping and thus lowers the effective mechanical frequency of the mirror j^. 
The effective frequency Ucs is given explicitly by 



2 2e^7Pir 



lu, + (7/2)2 

X ((7/2)^ -(^^-A^) .29) 

[(7/2)2 + (u;-A)2] [(7/2)2 + (^ + A)2]' ^ ) 

hence Ues scales as 1/ ^/l, that is, as the inverse square root of the mass. Note however 
that for the parameters of Fig. HI Ucs ~ O.Scu,^ for all choices of mass. Thus the maximum 
entanglement occurs for A/u^ ~ 0.8.) 

Physically, the condition A = uJcr corresponds to the optimum detuning for optome- 



chanical cooling 30(1 and can be understood in terms of cavity enhanced scattering of the 
anti-Stokes photons by the rotating mirror j4|. Thus it is not surprising that E]\[ reaches a 
maximum at this detuning. 



D. Angular momentum 

As pointed out in Ref. [l^, in the case of optical cooling of a vibrating mirror the 
linear momentum of the photons is difficult to change experimentally, and essentially entails 
replacement of the laser radiation source. For example to increase the linear momentum 
just by a factor of two would require an ultraviolet laser in place of an infrared laser. In 
contrast in the case of rotational cooling it is comparatively easier to achieve higher /. In 
light of this fact we also studied the entanglement as a function of the orbital angular 
momentum carried away from the mirror by the Laguerre- Gaussian mode,. The results of 
this analysis are summarized in Fig. [5l For the temperature T = lOK considered in the 
figure, entanglement appears only beyond the threshold value Ic ~ 21. This is consistent 
with Eq. Q which shows that for sufficiently low angular momentum, i.e. / < Ic, the 
optorotational coupling between radiation and the rotating mirror may not be high enough 
for entanglement to occur. The value of Ic depends on the other parameters that contribute 
to the optorotational coupling and of course also on the ambient temperature T. In general 
we found that entanglement appears when the effects of radiation torque coupling g are 
strong enough to counter the thermal noise, which increases perceptibly with temperature 
beyond T = T^. 
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FIG. 5: Logarithmic negativity Ei\f as a function of the orbital angular momentum carried by the 
Laguerre-Gaussian mode. The threshold angular momentum Ic ~ 21 at which entanglement begins 
to appear is indicated. Parameters of Table HI temperature T = lOK. 



No. 


Parameter 


Description 


Value 


Units 


1. 


L 


Cavity length 


1 


mm 


2. 


A 


Laser wavelength 


810 


nm 


3. 




Cavity resonance frequency 


27rl0i4 


Hz 


4. 




Rotating mirror angular frequency 


27rl0 


MHz 


5. 


M 


Rotating mirror mass 


100 


ng 


5. 


R 


Rotating mirror radius 


10 


/im 


7. 


Q(i> 


Mechanical quality factor 


2x10^ 




8. 


F 


Optical finesse 


5x10^ 




9. 


I 


Orbital angular momentum 


100 


h 


10. 


p. 
^ 111 


Input laser power 


50 


mW 


11. 


A/w0 


Normalized laser detuning 


1 





TABLE I: Definitions and approximate values of some of the parameters used in the text. 
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V. CONCLUSION 



To conclude, we have demonstrated theoretically that radiation torque can serve to en- 
tangle a Laguerre- Gaussian cavity field with a rotating spiral phase element. We have shown 
that the physical basis of entanglement is the exchange of angular momentum between the 
cavity mode and the rotating mirror. The mirror-field entanglement was investigated as a 
function of temperature and shown to be measurable at the temperatures routinely achieved 
by cryogenic apparatus. 

The entanglement was also shown to remain at its maximum value until a critical temper- 
ature at which the peak of the thermal spectrum occurred at the frequency corresponding 
to the rotor quanta, and beyond which the random excitation due to the thermal photons 
degraded the entanglement. The behavior of the entanglement with respect to cavity de- 
tuning was shown to be consistent with the mechanism of optomechanical cooling, i.e. the 
entanglement was maximum at a detuning known to correspond to optimal cooling of the 
mirror. Further, the fact that the angular momentum of the light can be changed relatively 
easily justified examining the entanglement with respect to changes in angular momentum. 
The main result here is that the mirror-radiation coupling only becomes strong enough to 
overcome the thermal noise and yield entanglement beyond a critical value of the angular 
momentum. 

The measurement of the entanglement investigated in this paper can be carried out in a 
manner analogous to that proposed in the case of a linearly vibrating cavity and illustrated in 
Fig. 3 of Ref . |9| . In the present case the mirrors are replaced by spiral phase elements, and 
a second optical cavity can be formed by placing a third transmissive spiral phase element 
beyond the rotating mirror of Fig. [TJ By using a weak intracavity field for this second cavity 
all the elements of the correlation matrix C can be obtained as proposed in Ref. and 
demonstrated in a measurement of the entanglement between two optical fields 
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APPENDIX A: LYAPUNOV EQUATION 



Here we derive Eq. (l24|) . which is a typical Lyapunov equation encountered in the study 
of dynamical systems. Also, we show our method for solving the equation. 



Derivation 



Beginning from Eq.(l23 



C 



dsM{s)DM{s 



(Al) 



where M(s) = e and the matrices B^D are given by Eqs.( fT3|) and ( l22l) respectively, we 
integrate by parts using /'(s) 



and gi^s) = De^^^ in the rule 



dsf\s)g{s) = [f{s)g{s)] 



oo 




dsf\s)g\s) 



(A2) 



This yields the equation 



C 



B-^e^'De^ 



dsB~'e''-We''^-'B^ 



poo 

-B^^D-B~^{ dse^'De 
Jo 

-B-^D - B-^CB'^, 



)B' 



(A3) 



where in the second step we have used the fact if the Routh-Hurwitz criterion is satisfied then 
lim e^' = 0, since all the eigenvalues of B have negative real parts. We then pre-multiply 

s— >oo 

the last line of Eq. ( 1A3I) by B and rearrange to obtain 

BC + CB^ = -D, (A4) 



which is Eq. ([2 



2. Solution 

Here we record our method of solving Eq. (1A4[) for the correlation matrix C. By definition 
[see above Eq. (!2Ti) ] C is a 4 x 4 real symmetric matrix. Thus we can write 

^ Al Oi a2 as \ 

ai X2 bi 62 

02 bi A3 Ci 

\ as 62 ci A4 / 



C 



(A5) 
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where we have denoted the diagonal elements with Greek letters and the off-diagonal with 
Roman letters. We now define a matrix C where 

C' = BC + CB^ + D, (A6) 

which equals zero as per Eq. fl2^ . Evaluating the RHS of Eq. flA6l) we find that at least one 
matrix element of C involves a single non-diagonal element of C. This observation, which is 
central to our method, turns out to be true at every iteration, where the process of iteration 
is defined below. 

We now solve for the non-diagonal element of C mentioned above by equating the relevant 
matrix element of C to zero. We then use that value of the non-diagonal element in C and 
re-evaluate C. This process is performed iteratively until all the off-diagonal elements of C 
are expressed in terms of the diagonal elements of C. At this stage C has the form 



/ \ n (Ai-A2)i^^ 7(A3+A4-1) \ 

^1 G 2G 



C 



„ , 7^(2n-2A2+l) 7^(A3+A4-l)-4AA2a;^-4Ai(G2-Aw^) 

(Ai-A2)u)^ 70(2n-2A2+l) > 7(2A3-1) 

G 2G 4A 

7(A3+A4-1) 7^(A3+A4-l)-4AA2aj0-4Ai(G^-Aa;^) ^(2A3-1) a 

V 2G 4Guj^ 4A ^4 

(AT) 

Also we find C = yields 4 equations in the diagonal elements Ai, 2,3,4 of C, and these can 
be easily solved, giving us analytical expressions for all the elements of the C matrix. 
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